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ABSTRACT 

In this paper we present the analysis of two broad band (0.1-150 keV) Bep- 
poSAX observations of the Seyfert 1 galaxy NGC 3516. The two observations 
were taken 4 months apart, on 1996 November 8 and 1997 March 12 respectively. 
We report a dramatic change in the degree of obscuration of the central source 
between the two observations, and propose, as possible explanations, transient 
absorption by either a stationary-state cloud of cold gas crossing the line of sight, 
or a varying-state, initially neutral and dense amount of expanding gas with de- 
creasing density and therefore decreasing opacity. 

We also report the detection of a second highly ionized absorber/emitter, 
which causes deep FeXVII-XXII K edges at ~ 7.8 keV to appear in both Bep- 
poSAX spectra of NGC 3516, and possibly produces the soft X-ray continuum 
emission and the 1 keV blend of Fe L recombination hues detected during the 
epoch of heavy nuclear obscuration. 

Subject headings: galaxies: Seyfert - galaxies: individual (NGC 3516) - X-ray: 
galaxies 

1. Introduction 

NGC 3516 is a bright (Fa-iofeeV = 7.8 x 10"^^ erg cm-^ s'^, Reynolds 1997) nearby (z 
= 0.009, D = 54 Mpc, for Hq = 50 km Mpc"^) Seyfert 1 galaxy which is known to host 
complex and variable systems of both mildly and highly ionized absorbers in the UV and in the 
X-rays. In the UV band, lUE spectra of NGC 3516 taken between 1978 June and 1989 October 
(Ulrich and Boisson 1983, Voit et al. 1987, Walter et al. 1990, Kolman ct al. 1993) show the 
presence of at least two systems of CIV resonant absorption lines: one is broad 2000 km s~^) 
and variable, the other is narrow (~ 100 km s~^) and constant. Both are produced by ionized 
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and outflowing gas which obscures both the continuum and at least part of the Broad Emission 
Lines (BEL). Four years later the broad absorption line component was no longer detected in the 
spectra of NGC 3516 taken with lUE (from 1993, Koratkar et al. 1996), or later with HUT ( in 
1995, Kriss et al. 1996a) and GHRS (in 1995, Crenshaw, Moran & Mushotzky 1998, Goad et al. 
1999, hereinafter GEA99), suggesting either a change in the covering factor or in the ionization 
structure of the gas responsible for that feature, or an increase in the emissivity of the CIV BEL. 
The complex X-ray absorber was observed by ROSAT (in 1992, Mathur, Wilkes & Aldcroft 1997, 
hereinafter MWA97) and ASCA (Kriss et al. 1996b, Reynolds 1997, George et al. 1998) through 
the detection of two strong absorption features at the energies of the OVII and OVIII K edges 
(0.74 and 0.87 keV respectively). The optical depths measured in the ROSAT and ASCA spectra 
using phenomenological 2-edge models, are similar to one another {tqvii ~ 0-6, toviii ~ 0-3). 

Neutral absorption has not been unambiguously found in this source. While the 1979 
Einstein observation (Kruper, Urry k. Canizares 1990) did not detect absorption by neutral gas, a 
subsequent low energy resolution EXOSAT observation in 1985 showed the presence of cold gas, 
with Nii^ ~10^^ cm~^ aS^ecting the low energy spectrum (Ghosh &; Soundararajaperumal 1991). 
A 1989 Ginga observation (simultaneous to lUE, Kolman et al 1993), showed a low energy cutoff 
of the intrinsic power law, indicating heavy absorption of the nuclear continuum along the line 
of sight. However, due to the limited energy coverage of Ginga (2-20 keV bandpass) it was not 
possible to establish the ionization state of that absorber, which was consistent with either a high 
column of neutral gas {N^"^'^= 4x10^^ cm~^) or a lower column of more highly ionized gas. 

In this paper we present data from two broad band BeppoSAX (Boella et al. 1997a) 
observations of NGC 3516. Preliminary results were presented by Stirpe et al. (1998). The paper 
is organized as follows. In §2 and §3 we present our analysis of the data, while in §4 we discuss 
our main results. Finally, we present our conclusions in §5. 
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2. Data Reduction and Analysis 

NGC 3516 was detected with a good signal-to-noise ratio in the three main Narrow Field 
Instruments (NFI) on BeppoSAX, spanning 0.1-120 keV. In Table 1 dates and exposure times of 
the two BeppoSAX observations of NGC 3516 are presented, along with the source count rates 
measured in the LEGS (Parmar et al. 1997), MEGS (Boella et al. 1997b) and PDS (Frontera 
et al. 1997). Data from these three instruments were screened following standard criteria, as 
detailed in Fiore, Guainazzi &; Grandi (1999). In particular, PDS data were screened using fixed 
rise-time thresholds. Data from the three MEGS units were merged to improve the statistics. 
The scientific products for the analysis are extracted from the equalized event files, using both 
the standard FTOOLS package (vs. 4.2) and the new release (vs. 1.1) of the ClAO (Chandra 
Interactive Analysis of Observation) software (Elvis et al. 2000). LEGS and MEGS source counts 
were extracted from circular regions of 6' radius, centered on the source. The same extraction 
regions were used to extract background counts from the long exposure "blank-fields" provided 
by the BeppoSAX Science Data Center^. PDS source and background counts were collected, 
respectively, during the on- and off-source modes of the instrument. For all the instruments, 
spectra were binned following the instrumental resolutions and allowing a maximum of 3 channels 
per resolution element. 

2.1. Model-Independent Evidence for Long-Term Spectral Variability 

NGG 3516 underwent moderate (factor of 1.5-2) 0.1-10 keV flux variations during both 
BeppoSAX observations (on a timescale of hours), however the ratios of the soft (0.1-2 keV) to the 
hard (2-10 keV) lightcurves were consistent with being constant, indicating no significant spectral 
variability across these energy bands on a 1-day timescale (the duration of each observation). 

A dramatically different behaviour was observed on the longer 4 months timescale between 
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the two observations, during which the source experienced strong spectral variabihty. The 0.1-2 
keV count rate increased by a factor of ~4, while the corresponding increase in the 2-10 keV 
band was only of 1.9 (Table 1). Figure 1 shows the model-independent ratio between the raw 
LEGS (0.1-2 keV), MEGS (1.5-10 keV) and PDS (13-100 keV) spectra for the two BeppoSAX 
observations (1996 November / 1997 March). The amplitude of non-linear effects in the responses 
of the NFIs (potentially affecting the results of such an analysis) is, at all energies, much smaller 
than the amplitude of the features visible in this ratio (Boella et al. 1997b, Parmar et al. 1997, 
Frontera et al. 1997). 

In the ratio in Figure 1, the 1-5 keV energy band is a smoothly increasing function of energy. 
It then flattens at E>5 keV, around a value of ~ 0.6. In the PDS band the ratio is ~0.8. This 
implies a steepening of the ratio between 10 and 13 keV. A sharp emission feature around 1 
keV is visible superimposed on these smooth changes. All this strongly suggests that some, and 
perhaps all, of the several components that make up the 0.1-100 keV spectrum of NOG 3516 
varied independently of each other from one observation to the other. In particular, the smooth 
and monotonic rising of the ratio up to ~ 5 keV, along with its flattening above this energy and 
up to 10 keV, suggest a drastic decrease in the degree of obscuration of the central X-ray source 
between the two observations rather than an intrinsic variation of the slope of the primary X-ray 
continuum. In this case the asymptotic value of ~ 0.6 around which the ratio flattens above 
~ 5 keV, would be the intensity increase factor of the intrinsic X-ray power law between the two 
observations. The higher PDS count ratio of ~0.8 may suggest a higher relative intensity of a 
Compton-reflection component in the 1997 observation. There is no clear indication for a similar 
effect at the energy of the Fe K line. The emission feature at 1 keV may be the signature of iron 
L emission, visible only during the 1996 observation, and therefore either external to the nuclear 
environment, and detectable by the BeppoSAX LEGS only when the soft nuclear continuum was 
heavily obscured, or internal to that, but only partially obscured by the amount of neutral (or 
mildly ionized) gas obscuring the X-ray source in 1996. Intrinsic variation of the intensity of the 
Fe-L complex in the hot plasma is, of course, a further possibility. 
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We use this scenario to guide the modehing and fitting of the spectra in the fohowing sections. 

3. Spectral Modelling and Fitting 

We have performed global modeling and fitting of the 0.1-150 keV spectra of both BeppoSAX 
observations of NGC 3516 using the XSPEC 10.0 package. In Table 2 we summarize the results 
of this analysis. Errors are quoted at 90% confidence level, for one interesting parameter (i.e. 

= 2.71, Avni 1976). The cosmology used is Ho=50 km s^^ Mpc~^ and qo=0.5. For all the 
absorption (neutral and/or ionized) models that we use in this paper we adopt solar abundances 
as estimated by Grevesse and Anders (1989) and Grevesse and Noels (1993). Photoionization 
models for fitting purposes include photoelectric absorption but not resonant absorption and gas 
emission, and have been built following Nicastro et al. (2000). Finally, we use pure Compton 
reflection models built by running iteratively the routine pexrav (Magdziarz & Zdziarski 1995) for 
the inclination angle of 30" (consistent with the best fit interval measured, see §3.1.3) and a range 
of values of the slope T of the incident continuum and the e-folding energy Ec of the high energy 
exponential cut-off. For fitting purposes we always combine the Compton reflection model with a 
cut-off power law, linking the two F and to the same values, and leaving the two normalizations 
free to vary independently. 

We first fitted the 1996 November and 1997 March spectra (hereinafter N96 and M97) with 
a model consisting of a power law plus Galactic absorption {Nh= 3.4 xlO^'^ cm~^, Dickey &; 
Lockman 1990). Both fits are very poor {x^/dofiN9Q) = lU7/94,x^ /do f{M 97) = 600/101), 
and the best fit photon indices are different: F ~ 1.2 and 1.5, respectively. The residuals clearly 
show the presence of several additional non-modelled spectral features (Figure 2a and b) . Some 
are common to both spectra: (a) a line-like feature at ~ 6.4 keV, (b) a deficit of counts at ~ 8 
keV, and (c) a systematic deviation across the PDS band. However, at low energy, (E<4 keV) 
the shapes of the residuals are quite different from one spectrum to the other. Our best fitting 
parameterizations of N96 and M97 are thus quite different, and are summarized in Table 2. 
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The two Best Fitting models (hereinafter 96BF and 97BF respectively) include several 
common components, which make up our "Base Model": (a) a power law with a high energy 
cut-off; (b) a neutral absorber (to account for Galactic absorption along the line of sight); (c) 
a hot photoionized absorber (HA, see Table 2), parameterized by a column density Nh and an 
ionization parameter U (to account for the 8 keV deficit) ; (d) a Gaussian emission line (to account 
for Ka iron fluorescence); (e) a pure Compton-reflection component (limited to F and Ec, as noted 
earlier). However, N96 and M97 differ from each other in requiring additional components that 
are not held in common: (f) a neutral absorber only partially obscuring the direct view of the 
primary X-ray power law and (g) soft emission by an optically thin plasma, for the N96 spectrum, 
and (h) a second, more mildly ionized absorber, for the M97 spectrum. 



3.1. The "Base Model" 

Our estimates for the parameter values of our "Base Model" were consistent with each other 
between the 96BF and the 97BF parameterizations (see Table 2) . In the following we describe and 
discuss the individual components of this model referring to the numerical results obtained for the 
M97 spectrum, the one with higher statistics. 



3.1.1. The Intrinsic Power Law and the High Energy Cut-Off 

The 0.1-100 keV primary continuum of NGC 3516 is well described by a simple power law 
F{E) = AE~'^ (F(E) in ph s'^ cm'^ keV"^) with T = 2.04+[|;[j^ (Table 2). This component varied 
in flux by ~ 70% between the two observations (i.e. A(N96)/A(M97) = 0.6), without changing in 
shape. We then conclude that little, if any, of the dramatic spectral variability experienced by the 
source can be due to variation of the primary continuum, either in shape or intensity. 

We could only estimate a lower limit (E> 450 keV) on the e-folding energy of the high energy 
cut-off, and were not able to measure any significant change of this parameter, between the two 
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observations. 

3.1.2. Iron Line and Compton Reflection 

Both the iron Hne and the Compton-reflection components were required by our data at high 
significance. Removing each of these components, one at time, from the 97BF model and refitting 
the data gives, respectively (with Au indicating the number of free parameters eliminated by the 
97BF model): Ax^ = -25 {Au = 3) and Ax^ = -41 {Au = 1). 

The best fit energy of the iron emission line (E = 6.41 it 0.15 keV) is consistent with that 
of Fel-XXII Ka transitions. When modelled with a simple Gaussian (as in our adopted best 
fitting parameterization 97BF), the line is consistent with being narrow, with a 90% upper limit 
on its width of 0.27 keV. We note that the energy resolution of the BeppoSAX-MECS is only 
of ~ 500 eV at 6.4 keV. Furthermore, an accurate determination of the exact shape of the iron 
line profile is hampered by the presence in the data of the ~ 8 keV edge-like absorption feature. 
However, we have checked the consistency of our data with more complex scenarios, replacing the 
Gaussian emission line in the 97BF model, with a relativistically broadened and distorted emission 
line from the accretion disk (the model diskline in XSPEC: Fabian et al. 1989), as suggested 
by previous ASCA observations (Nandra et al. 1997a, 1999). We fixed the internal radius on 
the accretion disk at 6 gravitational radii, the energy of the emission line centroid at 6.33 keV 
(the best fit energy in 97BF) and the spectral index of the emissivity law at -2, and refitted the 
data. The result was inconclusive. We obtained the same x^/i^ = 86/95 as in 97BF. The outer 
radius on the disk was only barely constrained and larger than 11.5 gravitational radii. All the 
other model parameters were consistent within the errors with the 97BF measurements. We also 
checked the extreme hypothesis that all of the negative residuals visible in Figure 2 at ~ 8 keV 
could be explained in terms of a very distorted blue-profile of the iron line, eliminating the HA 
component and refitting the data. This yielded a worse than 97BF (Ax^ = —18, for Au = 3), 
and accounted only partially for the absorption feature at ^ 8 keV. We then conclude that a 
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relativistically broadened and distorted emission line profile with no FeXVII-XXII K absorption is 
not a necessary component for the description of our data. 

Thanks to the broad band of BcppoSAX wc were able to deteet a bump in the 10-30 keV 
data (see Figure 2a and 2b). We model this bump as due to Compton reflection of the primary 
power law off a slab of almost neutral matter re-radiating the isotropic incident X-ray continuum. 
During both the 1996 and 1997 observations we measured relative fractions of reflection consistent 
with each other (Table 2), and fully consistent with the value of 1 corresponding to a solid angle 
of 0.5 as seen by a central isotropic source of primary radiation. Their ratio, 0.81q'2, is consistent 
with the value derived in §2.1 from our model independent analysis. If the same matter were 
also producing the Ka fluorescence iron line (but see also §2.1), we would expect the equivalent 
width of this line to be ~100-150 eV (e.g. Mushostzky, Done & Pounds 1993), consistent with the 
measured value (see Table 2). 

3.1.3. The Hot Absorber 

Both the N96 and M97 spectra of NGC 3516 show a clear absorption feature around 8 keV 
(Figure 2). A neutral reflection component alone cannot explain the edge feature at ~8 keV, even 
allowing the inclination disk parameter to vary in a wide range (from O" to 60"). We accounted 
for this feature by including a hot photoionized absorber (HA) in 96BF and 97BF, obscuring the 
line of sight. The same model was used for NGC 3516 to fit the absorption feature detected at the 
same energy (E~8 keV) by Ginga (Kolman et al. 1993). We stress again that the model used 
to fit the data includes neither gas emission, nor resonant absorption. Removing this component 
from both 96BF and 97BF and refitting the data, gives respectively Ax^ = —13 and Ax^ = —29, 
for Ai^ = 2. 

The best fitting parameters U and Nh do not vary (within the errors) between the two 
observations (Table 2), and their values indicate that a high column density {N^^ = (3.2^^^) x 10^^ 
cm~^, very similar to the value found by Kolman ~ 3 x 10^^ cm' -2) of highly ionized (U = 224tlf) 
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gas is absorbing the primary continuum of NGC 3516 in both data sets. In such an absorber 
elements hghter than oxygen and neon are fully ionized, and therefore do not imprint any visible 
features on the emerging spectrum. However, iron is still not fully ionized and contributes strongly 
to the opacity of the absorber at the energies of the photoelectric K and L absorption edges of 
FeXVll-XXII. Absorption by FeXVll-XXll K imprints deep absorption edges at ~ 7.8 - 8 keV on 
the emerging spectrum, while L absorption by the same species affects the whole 1 to 5 keV band, 
gently flattening a low-resolution spectrum rather than producing visible sharp and resolved single 
absorption features on it (unless Fc-L resonant absorption lines dominate on recombination lines 
^). Fitting low resolution data with a phenomenological single-edge model, accounting for the 
^ 8 keV absorption feature, would result in a wrong estimate of the slope of the intrinsic 0.1-100 
keV power law, which would appear to be flatter by AF = 0.28 it 0.13. Ionized, and possibly 
relativistically smeared, edge-like features at energies >7.1 keV have often been detected in the 
hard/low state of Galactic Black Hole Candidates, where the accretion disk is supposed to be 
hotter than that of AGNs by a factor ~30 (Done & Zycki 1999, Ross et al. 1996). In this contest, 
to further test the uniqueness of our interpretation, we applied an ionized reflection model, which 
possibly can also account for an ionized edge-feature, to our data. The Ax^ between a model 
with an ionized reflection (model pexriv in XSPEC, Magdziarz Sz Zdziarski 1995) and the best 
fit (BF) is Ax^/Az>=25/1 and Ax^/Azy=23/2 if the iron abundance parameter is left free to vary. 
The variation of the inclination angle of the ionized reflector {6 = 30°, 45", 60°) does not influence 
the test. Moreover, the expected ionized iron line at 6.7 keV, consistent with the ionization 
degree of the edge, is not detected, even though a ionized line could have been missed due to the 
low resolution instrument. 

The preferred interpretation is thus in terms of a high ionization, high column density 

^The relative contribution of Fe-L resonant absorption lines, compared to emission lines, depends 
critically on the geometry and dynamics: if the gas has a high microturbulence velocity and a lower 
or comparable bulk-velocity, then Fe-L resonant absorption dominates and imprints a visible "1- 
keV" feature on low resolution spectra (Nicastro, Fiore & Matt, 1999). See also §4.2 
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absorber which can simultaneously account for both the ~8 keV edge and the overall continuum 
shape at medium energies and, most important, it provides the statistically better description of 
the available data. 

3.2. Departure from the "Base Model" 

As pointed out in §3 the "Base Model" does not provide a fully satisfactory description 
of either the N96 or the M97 spectra. Additional components are present in our best fitting 
parametrizations 96BF and 97BF. We describe them in the following two sections. 

3.2.1. The "Cold" Absorber and the Soft Thermal Emission in N96 

To model the 0.1-5 keV residuals of Figure 2b, we added to our "Base Model" a neutral 
absorber of column density N^"^'^ which only covers a fraction Cgou of the nuclear X-ray 
source. Fitting the N96 data with this model produces a significative improvement in the 
{Ax^/Aiy = 73/2) and gives Ng"'*^ = 2.3l1:oil x 10^^ cm'^ and 0.7 < Csou < 1- (We note that 
using an ionized absorber model the ionization parameter is consistent to be zero, while the 
column density of the warm gas is ~ iV^°''^, indicating absorption by neutral gas as the most likely 
explanation). Nevertheless, the residuals continue to show a clear line-like unresolved feature 
around 1 keV (Figure 3, upper panel). 

At this energy, recombination L lines from FeXVll-XXIV are expected to strongly contribute to 
the emission from either a coUisionally ionized plasma (CIP) with a temperature of ^ 1 — 2 keV, 
or from a photoionized plasma (PIP) with a high ionization degree (U^ 200) and column density 
{Nh ~ 10^3 cm-2). Only for fitting purposes, we use here a CIP model, while a more accurate 
discussion on the physics of this component and on its possible identification is deferred to §4.2. 
We added a CIP component (the model mekal in XSPEC, Mewe, Gronenschild & van den Oord 
1985) to our model, and refitted the data. The fit is again significatively improved by the addition 
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of this component (Ax^ = 13, for 2 additional parameters: model 96BF), and the residuals are 
now flat over the entire 0.1-100 keV band. The best fit temperature of this collisional plasma 

is kT = 1.78^Q3(!) kcV. The covering factor of the neutral gas is now larger (Csou > 0.95) than 
in the previous fit, indicating that only a small or null fraction of the nuclear radiation escapes 
absorption. 

3.2.2. The "Warm" Absorber in M97 

Our "Base Model" did not provide a satisfactory description for the M97 spectrum either. 
Several features remained at E<2 keV as shown by the residuals in Figure 3b (bottom panel). 
These strongly suggest the presence of a typical "warm" absorber obscuring the line of sight and 
imprinting deep OVII and OVIII K edges, at 0.74 and 0.87 keV (see Reynolds 1997, George et al. 
1998). For this source the signature of a complex "warm" absorber had been previously detected 
by ROSAT-PSPC and ASCA (MWA97, Kriss et al., 1996b). We then added a photoionized 
absorber component to our "Base Model" and refitted the M97 data. The residuals are now flat 
over the entire BeppoSAX band, and the reduced is 0.90 (Table 2: 97BF model). 

4. Discussion 
4.1. The Vciriable Absorber 

The 1996 November spectrum of NGC 3516 is heavily absorbed by almost neutral gas which 
covers most of the primary continuum source, (here we choose 

Csou ~ 0.95, the best fit value 
obtained fitting the soft emission with the CIP component, §3.2.1), and whose equivalent hydrogen 
column density (as derived by a fit with a neutral absorber) is as large as N^"'''^ = 2.3 x 10^2 cm 
(exceeding the Galactic value by almost two orders of magnitude). This value would be even 
larger if the gas were mildly ionized (OII-III dominating). There is no evidence for variation of 
the degree of obscuration of this gas during the observation (~20 hours) and this component is 
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not required by the 1997 March data (4 months later) ^, and the UV continuum of this source in a 
HST/FOS observation taken ~ 18 day after the 1996 November BeppoSAX observation (GEA99) 
shows no marked extinction by dust (assuming a Galactic gas/dust ratio we should expect A^, ~ 10 
mag, for Nh = 2.3 X 10^2 Gorenstein, 1975). 

All this allows us to put constraints on the geometry and the physical state of this component, 
and to speculate on different possible scenarios. 

4.1.1. Transient Cold Absorption by a Broad Emission Line Cloud? 

The typical timescale of variability of the X-ray continuum (e.g. Nandra et al. 1997b) gives 
an upper limit on the linear size of the emitting region: D ~ cAt ~ 10^^ cm. In order for an 
homogeneous cloud of absorbing gas to cover almost the entire primary emitting region, either the 
cloud is at a distance R close to the central source, comparable to D, and its linear sizes I are 
comparable and/or larger than D, or R is much larger than D and the source appears point-like 
as seen by the cloud. Photoionization arguments strongly support the second picture. In fact, 
assuming R ^ D gives an ionization parameter U> 2500 , where nio is the electron density of 
the gas in units of 10^° cm~^. We used the ionizing luminosity of NGC 3516, derived assuming 
the 2200 A flux and aox from George et al. (1998) and X-ray fluxes and spectral shape from our 
data. Such an ionization parameter would be associated to almost neutral gas (U^O.Ol: oxygen 
distributed between OI-OIII) only for implausibly high electron densities ^ nio > 2.5 x 10^. More 
plausible electron densities of nio = 0.1 — 100, as those estimated for the Broad Emission Line 
Clouds (BELCs) in AGNs (Netzer 1991, Ferland et al. 1992, Krolik et al. 1991), would result in 

^We note, however, that an ionized, similar Nh, absorber is needed in 1997: see §4.1.2. 

^We note that this is one of the highest densities one may expect to find in the nuclear environ- 
ment of a Seyfert 1, being of the order of the density of the matter in a standard Shakura-Sunyaev 
(1973) accretion disk around a black hole of 10'' — 10^ M©, accreting at 0.1-0.01 times the critical 
rate, and in proximity of the last stable orbit. 
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very highly ionized gas, which is inconsistent with our data. We then assume that R» D. 

Reverberation mapping studies of the BELs in NGC 3516 (Wanders et al. 1993) have 
estimated a distance Rbelc of the BELCs from the source of ionizing continuum, of ^ 11 
Zc/ (~ 3 X 10^^ cm), supporting the possibihty that the gas obscuring the X-ray continuum of 
NGC 3516 in 1996 is associated with a BELC. At this distance, and assuming nf^f ^'^ = 0.1 - 100, 

the ionization parameter measures V{Rbelc') ^ 0.0003 — 0.03, and the gas is only mildly ionized 
and therefore compatible with the nature of the absorption seen in the X-ray band in 1996. 

All this suggests that a BELC crossing our line of sight could have caused the obscuration of 
the X-ray primary source in 1996. As a consistency check we can calculate the time tcross needed 
for a BELC to cross our line of sight. Assuming as transverse velocity vbelc, the measured H/? 
FWHM of 4,000 km s~^ (Wanders, et al, 1993), we obtain tcross = l/v ^ 2, 500iV^coicinj;o^t;4000 
(where I is the linear size of the cloud). For the measured N^°''', we have then tcross ~ 6, OOOnJ^Q^ 
s. Since we know that tcross must be longer than 20 hours (the duration of the 1996 observation) 
and shorter than ^ 4 months (the time elapsed between the two BeppoSAX observations), the 
above equation constrains the electron density of the gas to the range 6 x 10"*^ cm~^<raio^ 0.1, 
We stress again that, if the gas were mildly ionized (as is the case for the BELCs), then we would 
be underestimating the equivalent N^"'*^ obtained fitting the data with a neutral absorber model 
and therefore the interval of densities given above would shift towards higher values. This range of 
densities contains the estimated BELC value. A much narrower interval could in principle be set 
using the time elapsed between the 1996 November observation and the last HST/FOS observation 
of GEA99), showing no Mgll resonant absorption nor extinction in the UV continuum, in At = 18 
days. This would give a lower limit on the density of nio^O.Ol. 

In the framework of the proposed scenario, a BELC crossing the line of sight, we estimated 
the probability of such an event, based on an estimate of the filling factor of the BELCs. The 
total volume of a sphere with radius Rbelc is Vtot ~ 8.2 x 10^^ cm^. From the luminosity of the 
broad H/3 of NGC 3516, l^np ~ 2.3 x IQ'^^erg s'^ (Rosenblatt et al. 1994), we can estimate the 
effective volume occupied by all the BELCs, given by: 
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Veff = LHp/(x{Te)nl. For an electronic temperature Tg of ~ 10^ K, at the energy of the H/3 
emission line, the recombination coefHcent a{Te) is ~3.4xl0~^^ erg cm^ s~^ (Osterbrock 1989) and 
Ue ~ W^^cm~^ (Netzer 1991). We then obtain K;// ~ 6.7 x 10^"^ cm^. The probability for a cloud to 
find itself along a particular line of sight at any radius is then given by P ~ Veff/Vtot ~ 8.1 x 10~^. 
This low probability supports the picture of an event which occurs very rarely, as indeed observed 
in AGN (Kruper, Urry & Canizares 1990). 

4.1.2. A Variable State Absorber? 

Another intriguing possibility is that the ionization state of the absorber varied between the 
two BeppoSAX observations, becoming more ionized, and therefore transparent to the soft X-ray 
radiation. In this scenario we may speculate that the "warm" absorber seen in the 1997 March 
spectrum of NGC 3516 is the same gas which was obscuring the line of sight 4 months earlier, but 
with a very different ionization degree. The same picture was drawn for this source by MWA97 
who proposed a model for a variable "warm" absorber in NGC 3516 to reconcile the UV and the 
X-ray spectra within a unique scenario (see also next section) . In the MWA97 model the presence 
and subsequent disappearance of a broad CIV absorption component at 1549 Awas consistent 
with a variable state X-ray absorber, from mildly ionized during the 1989 Ginga observation to 
highly ionized in subsequent ROSAT-PSPC and ASCA observations. 

In the framework of the scenario above described, we derived some dynamical properties 
of this variable state absorber. At a given ionizing flux, the ionization degree of a photoionized 
absorber can vary due to variations of the distance of the gas from the ionizing source, or of the 
electron density in the gas (or both). In the first of these two extreme hypotheses we assume that 
the absorber has a constant density in the range nio = 10-^ - 10-1 (typical of "warm" absorbers 
in Seyfert Is: Nicastro et al. 1999, 2000, MWA97). The gas is almost neutral during the 1996 
November observation, which gives an upper limit on its ionization parameter of U ^0.01, above 
this value, and for the ionizing continuum shape of NGC 3516, the opacity of the gas at low 
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X-ray energies drops drastically. This gives a range of possible distances from the central source 

cm (the larger the density the smaller the lower limit on R) . If the gas 
were moving toward the central source with a radial velocity Vrad, keeping its density constant, 
it should have required less than t<4 months to cover the distance AR = R/10 needed for its 
ionization parameter to increase by a factor of ^ 100 (to reach the values observed during the 1997 
November observation). The estimated radial velocity of the gas is thus Vrad ~ Ai2/t> 0.06c. This 
limit would become even higher if the cloud became denser while collapsing toward the central 
source. We note that the most direct implication of this extreme hypothesis is that UV and high 
resolution X-ray spectra of Sy 1 galaxies, during such relatively shorts events, should show highly 
redshifted resonant absorption lines with varying strenght. Such phenomenon has never been 
observed so far. 

Alternatively, the change in the ionization state of the absorber might be a result of a 
change in density (as in MWA97). To change from an almost neutral state (U^O.Ol) to a warm 
state {U ~ 5.37) would require a drop in density (Uoc njj^) by a factor of i;500. For example, 
the density of a BELC would have to change from a few 10^*^ cm~^ to about 2 x 10^ cm~^. A 
column density of 10^^ cm~^ and a density of 10^'^ cm~^ imply a thickness of 10^^ cm for the 
cold absorber. Similarly, the thickness of the warm cloud would be about 5 x 10^^ cm. Such an 
expansion has to occur on a time scale of $ 4 months. If the cloud is expanding adiabatically, 
its velocity dispersion would then have to be about 500 km s^^. Such a velocity dispersion is 
observed in associated absorption lines in AGN, so the derived value is completely reasonable. 
We conclude that a scenario based on a variable state absorber is plausible (see §4.3 for further 
details). Unfortunately, no other X-ray observation was performed between the two BeppoSAX 
pointings to further constrain the variability time scale in the soft X-ray energy band. 

Both cases discussed above were based on the assumption of no changes of the intensity of the 
photoionizing continuum. Variations of the ionizing flux could have contributed to the observed 
changes of the ionization degree of the gas. If, between the two observations, the continuum 
increased by a factor of ~ 100 and then recovered the level measured during the 1997 observation. 
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then the gas could have been almost instantaneously ionized by this dramatic event, but not 
yet had enough time to recombine to its equilibrium neutral state. Alternatively, the flux level 
seen during the 1996 observation could have been preceeded by a long duration event of very low 
intensity which caused the gas to recombine to an almost neutral state. The gas during the 1996 
observation would then be underionized with respect to the observed ionizing continuum level. 
Both cases would require low electron densities of the absorber (rie^lO^ cm"^, see Nicastro et al. 
1999, for details). 



4.2. A hot Ionized Absorber/Emitter gas? 

Two apparently unrelated spectral features are clearly visible in the N96 spectrum of 
NGC 3516: (a) soft emission (with a hump at ~ 1 keV), not obscured (or only partly obscured) 
by the amount of neutral gas covering the primary X-ray source, and (b) the signature of heavy 
absorption by highly ionized gas at ~ 8 keV. In our spectral analysis of N96 we modeled these 
two features as due to two distinct components: a collisionally ionized emitting plasma and 
a photoionized absorber, respectively. Extended emission by the host galaxy surrounding the 
AGN has frequently been observed in Seyferts: e.g. Elvis ct al. (1990), Wilson ct al. (1992), 
Weaver et al. (1995). However, in the particular case of NGC 3516, the extended component 
was probably unresolved by the ROSAT-HRI (Morse et al. 1995). Wilson et al. (1996) found 
that in their sample (which includes NGC 3516) the 0.1-2.5 keV intrinsic luminosity of the 
extended components ranges between lO^'^ and 10^^ er^ s~^. In the present data we measure 
Lo.i-2.5keV = 5 X 10^^ erg s^^, a value only marginally consistent with the observed range, even 
allowing for intrinsic variability. Based on the correlations found between far infrared (FIR) and 
2 keV luminosity of bright spiral galaxies (Fabbiano et al. 1988), we also compared the 2 keV 
luminosity L2keV = 9 x 10^-*^ erg keV""^ measured, with the mean L2keV value for the 11 
objects of the Fabbiano et al. (1988) sample with Lpm = 10^^-^ — 10^^-^ erg s~^. This range 
includes the value of 10*^^ erg obtained for NGC 3516 (Bonatto & Pastoriza 1997). We found 
^2keV — ^^^^ ~ W^^erg s~^, again much lower than measured in our data. We then conclude that 
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the soft X-ray emission visible in 1996, cannot be entirely due to an extended component in the 
host galaxy. 

An intriguing alternative is that both the ~ 8 keV absorption edges visible in N96 and M97 
and the soft emission in N96 could be self-consistently accounted for by a unique photoionized 
absorber /emitter in the nuclear environment of NGC 3516. If photoionized, this gas has a 
temperature of T ~ 2 x 10® K, and iron (the lightest abundant element not yet fully ionized) is 
smoothly distributed between FeXVll and FeXXVI. Due to the high temperature the emissivity of 
this gas is large and contributes to the overall spectrum in an amount which depends linearly on 
its covering factor around the central source. In particular, a strong contribution by FeXVII-XXIV 
L recombination lines is expected at ~ 1 keV, where a clear emission feature is present in the N96 
spectrum. Models for photoionized gas accounting for gas emission, other than photoelectric and 
resonant absorption, can then be used to constrain the geometrical configuration of the observed 
absorber / emitter. 

In the following wc try to test this hypothesis quantitatively. The basic assumptions of this 
hypothesis are: (a) the neutral absorber in 1996 covered almost entirely {Cgou > 95 %, as from 
the fit with the CIP component, §3.2.1), the nuclear source of primary X-rays but left uncovered 
a much larger fraction (1 — Cenv) of the extended nuclear environment (i.e. the warm and hot 
scattering media) which is seen through scattering of the primary continuum; (b) during the 1997 
observation both the primary nuclear continuum and the reprocessed components are seen directly 
(c) no other component can significantly contribute to the line emission. 

We proceeded to build (following Nicastro et al. 2000) two series of photoionization models, 
which include gas emission and resonant absorption, for the two ionized absorbers observed in 
NGC 3516. We used column densities and ionization parameters as measured fitting the N96 
and M97 spectra with models accounting for photoelectric absorption only (Table 2). We built 
these models for a number of values of the covering factors fwarm and fhot as seen by the central 
source, varying from to 1, and for particular values of outflowing (vout) and microturbulence (a) 
velocities of the absorbers/emitters. In particular we used t'^arm = ^00 km s"'^, (Twarm = 1; 100 
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km s-^ (observed in the UV for the broad CIV system of NGC 3516, Kolman et al. 1993), 
and Vf^^ = 500 km s~^,ahot = 200 km s~^. For the "hot" absorber the adopted ratio v^^l/ahot 
maximizes the contribution of line emission compared to resonant absorption (Nicastro et al. 
2000). The particular choice of the absolute values of v and a, separately, is however not critical. 
The warm component was not included in the models for the N96 spectrum, since we did not 
see it in the data (which can be either because the component was not present or because it 
was completely obscured by the amount of almost neutral gas covering the line of sight during 
that observation). We also neglected the emission contribution from this gas. However given the 
temperature and ionization degree measured in 1997 (in the photoionization hypothesis), emission 
of the warm component is expected to be very weak, and should influence mostly the 0.5-0.6 keV 
portion of the spectra where weak OVII and OVIII emission lines are expected. This contribution 
would be hardly detectable with BeppoSAX (Nicastro et al. 2000). The hot component has been 
instead included in both models. We held fhot constant between the two observations. For the 
"cold" absorber in 1996, we used the best fit column density and Cgou from Table 2. 

We did not fit these models directly to the N96 and M97 spectra, but proceeded as follows. 
We varied the values of the three parameters fhot, fwarm and Cenv to let (a) the 0.1-2 keV flux 
predicted by models for N96 match the best fitting 96BF 0.1-2 keV flux, and (b) the 0.1-10 
keV ratio between N96 and M97 models match the observed raw data ratio N96/M97. In this 
particular scenario, we could put relatively strong constraints on the covering factors fhot and 
Cenv, but not on fwarm, which we then fixed to the value of 0.5. We found that fhot and Cenv 
lie in the ranges 0.3-0.4 and 0.15-0.25 respectively. Figure 4 and 5 show two models for the N96 
and M97 spectra of NGC 3516 (Fig. 4a,b, upper and bottom panels respectively), and their ratio 
(Fig. 5) in the energy band 0.1-10 keV, for fhot = 0.35 and Cenv = 0-2. The ratio between the raw 
0.1-10 keV N96 and M97 data is also plotted in Figure 5c. The agreement between the models 
and the data is quite good. In particular we note that the 1 keV emission feature present in the 
data is well reproduced by Fc L emission by the hot photoionized components, which was much 
more visible in 1996 due to the heavy obscuration of the primary X-rays by the "cold" absorber. 
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We then conclude that a highly ionized component is present in the nuclear environment of 
NGC 3516, and that its emission is visible only when our view of the primary X-rays is at least 
partially covered by a large amount of neutral absorber. 

This component could well be present in all Seyferts, but has eluded detection so far. If this 
component is spherically distributed around the central source, and if the covering factor fhot and 
column density N^°'' measured in NGC 3516 arc representative of the entire class, then we would 
expect that absorption K edges due to highly ionized iron should be observed in almost 30-40% of 
all known Seyferts. Furthermore, in all heavily obscured sources (i.e. Seyfert 2s, or objects with 
transient absorption, like NGC 3516), we should detect emission by this hot gas, at energies below 
1-2 keV. Soft X-ray emission is indeed frequently observed in many bright and well studied Scyfcrt 
2s (Turner et al. 1997, Comastri ct al. 1998). High resolution spectrometry, like that available 
with Chandra and/or XMM, will test this speculation unambiguously by allowing the detection of 
narrow emission and/or absorption lines even in objects with a low degree of nuclear obscuration 
with high accuracy, as already studied in NGC 5548 and NGC 3783 (Kaastra et al. 2000, Kaspi et 
al. 2000). 

4.3. The Complex X/UV Connection 

Following MWA97, we investigate whether the X-ray absorbers in NGC3516 will produce any 
associated UV absorption lines. 

If the cold absorber is mildly ionized, the maximum degree of ionization would be ?7 ~ 0.01; 
as seen in §4.1.2, a higher value, with the ionizing continuum shape of NGC 3516, would imply a 
drastic drop of opacity of the gas at low X-ray energies. In such an absorber, the dominant ions 
would be OI-III, Mgll-III and CII-IV. For a total column density of a few times 10^^ cm~^, the 
predicted column densities of Mgll and CIV are 2.5 x 10^^ and 4.7 x 10^^ cm~^ respectively. For 
the warm absorber, the CIV column density would be 5.4 x 10^^, but Mg is too highly ionized to 
have observable Mgll. For the cold absorber, the Mgll and CIV column densities are very large 
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and even for a velocity dispersion parameter (b) of 100 km s ^, the resonance absorption lines 
would have large equivalent widths (EW ^ SA) . 

In the warm phase, the CIV column density is less than that in the cold absorber. However, 
the predicted EW of the CIV A1549 absorption line would still be detectable (EW>1.0 A, for 
b=100 km s~^). The absorption line EWs would be even larger for larger b. Therefore, for 
reasonable parameters, the X-ray absorbers in NGC 3516 will contribute to the absorption in UV. 
How broad the lines would be will again depend upon the b parameter. The UV observations 
closest in time to our X-ray observations {At = 18 days) are presented by GEA99 (§4.1). With the 
HST FOS, intrinsic Mgll absorption was not detected. This constrains the cold X-ray absorber (if 
still along the line of sight during the UV observation) to be almost completely neutral. GEA99 
found CIV A1549, 1551 absorption doublets with FWHM 651±16 km and 405±9 km s"^ 
respectively and EWs 2.87±0.09 A and 1.47=b0.04 A respectively. The lines are not resolved with 
the FOS and are likely to be saturated. GEA99 quote a CIV column density in the range between 
6.3 X 10^^ and 10^^ cm~^, very similar to our predicted value for the warm absorber. Earlier 
observations with GHRS have resolved the CIV absorption line into four distinct components, 
two of which are from the host galaxy of NGC 3516. The two nuclear components have column 
densities consistent with the GEA99 values. We conclude that the nuclear CIV absorption line 
(possibly component 1 of Crenshaw et al. (1998)) is likely to arise in the warm absorber observed 
by BeppoSAX. At the least, the warm absorber makes a substantial contribution to the absorption 
seen in the UV (based on CIV strength). 

The hot X-ray absorber is too ionized to contribute to the UV absorption lines, since carbon, 
oxygen and magnesium are fully ionized. 

X-ray and UV observations show that the system of X-ray absorbers in NGC 3516 is clearly 
complex and highly structured with multiple ionization stages just as in the UV absorption lines 
(§1). This complexity is enhanced by variability on many different time scales, from days to years. 
It is possible that on the very short time scales of a day, the variability of absorption is governed 



-22- 



by the photoionization-recombination timescales, while on the longer time scales from months 
to years the dynamical time may play the dominant role, either via bulk-motion or adiabatic 
expansion. It is conceivable that the different ionization stages trace different phases in the 
evolution of an absorber. A neutral absorber may expand as it outflows, becomes a warm absorber 
and eventually becomes hot and later completely transparent in the X-rays, or the "components" 
may be part of a quasi-static structure in an outflowing wind (Elvis 2000, Murray &; Chiang 1995). 
Again high resolution Chandra and XMM observations of this source will be crucial to disentangle 
the many X-ray components, and clearly establish (or reject) a link with the UV components. 

5. Conclusions 

We presented two BeppoSAX observations of NGC 3516, performed in 1996 November, and 
1997 March. 

Our main findings are: 

• The source underwent strong spectral variability between the two BeppoSAX observations, 
taken 4 months apart, due entirely to the drastically different degree of obscuration of the 
central source. During the first observation the nuclear X-ray continuum was absorbed by 
an equivalent column of cold hydrogen of 2.3 x 10^^ cm~^, while in 1997 November the 
absorption was fully compatible with that from our Galaxy. This is the first time that such 
a change in the degree of obscuration of a Seyfert 1 has been clearly observed. 

We investigated several possibilities and propose as possible explanations either a BELC 
crossing the line of sight, or a state-varying absorber. 

• We discovered the presence of a very highly ionized absorber/emitter in the nuclear 
environment of NGC 3516. This gas is visible in absorption (and with a stable physical 
state) in both BeppoSAX observations. 

Thanks to the large fraction of primary X-rays absorbed by neutral gas in 1996 we also 
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were able to detect emission from this gas, and to give an estimate of its covering factor 
fhot = 0.3 - 0.4 
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Table 1. Observation log of the two BeppoSAX pointings of NGC 3516. Count rates are 
given in the bands 0.1-2 keV, 2-10 keV and 15-150 keV, for the LEGS, MEGS and PDS 

respectively. 



Dates 


Exposures 




Rate (cts/sec) 






LEGS MEGS PDS 


LEGS 


MEGS 


PDS 



8-11-1996 16059 56013 20758 0.034±0.001 0.36±0.003 0.72±0.04 
12-3-1997 16909 57712 21713 0.14±0.004 0.69±0.004 0.86±0.04 



Table 2. Best fitting parameterizations 96BF and 97BF. Both include: a cut-off power 
law plus Galactic absorption a gaussian line plus cold reflection and a highly ionized 
absorber. 96BF includes also a partial covering absorber {Cgou ~0.95), and a thermal 
component. 97BF includes instead a "warm absorber". See §3 for details. 
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Fig. 1.— Ratio of the N96 and M97 spectra of NGC 3516. 

Fig. 2. — Residuals of a simple power law plus Galactic absorption fit for N96 (upper panel) 
and M97 (lower panel). 

Fig. 3. — Upper panel: N96 residuals of the "base" model plus a partial covering neutral 
absorber; an emission feature at 1 keV is clearly present. Lower panel: M97 residuals of the 
"base" model; the typical warm absorber signature is visible. 

Fig. 4. — 96BF model with no CIP component (upper panel) and 97BF model (lower panel). 
Photoionization models include gas emission and resonant absorption other than photoelec- 
tric absorption (see text for details) . 

Fig. 5. — Ratios between the two models in the upper and lower panels, and the raw 0.1-10 
keV N96 and M97 spectra. 
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